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NOMENCLATURE* 
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Source strength coefficient 
Speed of sound 
Energy flow per unit area 
Mean energy flow per \init area 

•fl 

Wave number in x-dlrection (Cartesian coordinates) 

Wave number in y-directlon (Cartesian coordinates) 

Wave number = (Cartesian coordinates); 

Wave number in axial direction (cylindrical coordinates) 

Reduced wave number in x-dlrection = k./<o* (Cartesian 
coordinates) 

Reduced wave number in y-directlon = k^/ <u' (Cartesiaji 
coordinates) 

Reduced wave number = k/cj* = (Cartesiaji coordinates); 

Reduced wave number in axial, direction = k/<u* (cylindrical 
coordinates) 

Mach number 

an integer 

Pressure Increment from acoustic waves 
Mean-square pressure from acoustic waves 
Radial cylindrical coordinate 


*The secondary symbols not defined in this list are defined 
where they are Introduced in the analysis. 
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Xq Radius of uniform velocity jet core 

r^ Radius of jet 

r^ Radial position of source 

R Radial spherical coordinate 

Strouhal number 
t Time 

X Streamwise coordinate (coordinates fixed in source) 

Xg Streamwise coordinate (coordinates fixed in ambient air) 

y Lateral coordinate 

z Normal coordinate 

R.P, Denotes ”real part of” 

I.P. Denotes "imaginary part of" 

7] Radial displacement of a fluid particle 

0 Angular position in far field measured from upstream 

jet axis 

P Fluid density 

^ Velocity potential 

Angle about x-axis (cylindrical coordinates) 

CO Generating frequency 

CO* co / c 

CO* Observed frequency 

coq Transient source frequency 
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Nomenclature , continued 
Subscripts 

c refers to critical 

j refers to jet 

s refers to source 

R refers to real part 

I refers to Imaginary part 



THEORETICAL STUDY OF REFRACTION EFFECTS 


ON NOISE PRODUCED BY TURBULENT JETS 

E. W. Graham and B, B, Graham 
Graham Associates, Shaw Island, Washington 


INTRODUCTION 

Ohjectives 

1 

This report is an extension of the work published in NASA CR-2390 . 
None of this work is intended as a comprehensive theory for the predic- 
tion of noise produced by turbulence in jets. Instead we concentrate on 
one aspect of this complex problem, the transmission of acoustical dis- 
turbances from the interior of the jet into the ambient air. These 
acoustical disturbances are generated by mathematically defined point 
sources drifting with the local fluid. We neglect temporarily the more 
difficult problem of identifying these mathematical sources with the 
turbulence which they are intended to represent. 

Because of the diversity of the material presented here (the far 
field of supersonic jets as contrasted with the near field of subsonic 
jets) this report is divided into two parts. 

In Part 1 we study the peaks which arise in plots of mean-square 
pressure versus angle in the far field for some supersonic jet examples. 
Such peaks may occur even when the basic wavelength produced by the 
pulsating source is several times the jet thickness. These peaks axe so 
severe in some instances that a pressxire change of several orders of 
magnitude may appear and disappear within a one degree change of angle 
in the fax field. 

In Part 2 the near field disturbances created by transient sources 
in a subsonic jet axe studied. This is a preliminary attempt to get 



theoretical data for compaxison with the escperimental in-house research 
work of NASA TN D -7269 ty Maestrello. It may also he useful for check- 
ing the applicability of infinite- jet analysis to realistic expanding 
jets of finite length. 

Related Work 

This report and the preceding one, NASA CR-2390, utilize the radia- 
tion fields of sources convected with a moving fluid. Prior work of this 

3 ^ 

type was done by Gottlieb , and by Moretti and Slutsky . More recently 
an analysis of this type has been made by Mani'^, 

Work on diverging jets (which are of course more realistically 

shaped than the cylindrical jets of the above literature) has been done 

678 

by Schubert , by Liu and Maestrello' , and by Padula and Liu . We hope 

such analyses can be extended to cover sources convected with the moving 

fluid. 

Possibly we have overlooked other closely related work which should 

be cited here, and certainly we have failed to refer to many jet noise 

studies which bear some relation to our own work. No disparagement of 

such omitted references is intended. 

In our preceding report (NASA CR-2390) we should also have cited 

5 9 

theoretical analyses by Mani and by Morfey and Tanna , and experimental 

10 

work by Ingard and Singhal , when mentioning the existence of more than 
one type of moving source. 


Description of Model 

The noise generator chosen is a sequence of transient sources 
drifting with the local fluid. The jet (either two-dimensional or 
circulax cylindrical) extends to infinity upstream and downstream, with 
velocity profile independent of streamwise position. Thus the large 
velocity gradients across the jet axe accounted for, and the smaller 
gradients in the streamwise direction axe neglected. It seems reason- 
able to suppose that the major refraction effects will be shown by such 
a model. However investigation of the streamwise extent of near-field 
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disttirbances may also cast light on the suitability of this model, 

In making this analysis we consider that turbulence ("self -noise”) 
is the only true originator of noise, and that "shear noise", being 
composed of linear terms, is part of the transmission process. 

It must be emphasized that this is not a stability analysis. We 
deal with a distribution of turbulence in the jet which is essentially 
independent of time. This steady-state situation is the end result of 
the action of instabilities. 

We then attempt to find the effect of the jet mean velocity profile 
on the transmission of acoustic disturbances from one element of turbu- 
lence through the jet and into the ambient air. The scattering effect 
of other elements of turbulence is neglected. 

Description of Methods* 

y -y - 

GeneraJ-ity of Methods 

Nearly all of the work in this report and the preceding report, 
NASA CR- 2390 , has in common certain necessary steps and concepts. The 
noise may be generated by sources (i.e, monopole sources), dipoles or 
quadrupoles. The noise generators may be on center, off center in a 
uniform velocity region or, in some cases, in the shear layer. Temper- 
atures may be ambient or varying across the jet. The jet itself may be 
two-dimensional or circulax. Still, with only minor modifications cer- 
tain basic ideas apply, and we review them briefly. 

For convenience we will describe the jet as being two-dimensional. 
The velocity profile must be independent of streamwise and lateral 
positions. For simplicity we will say "source", and speak of a single 
(constant velocity gradient) shear layer on each side of it, but other 
singularities can be used, and more complicated shear layers can be 
treated. 


*Much of this is reproduced from NASA CR-2390 for the reader's conven- 
ience. 

**See references 11-15. 
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The Source in a Jet 


It is first assumed that the pulsatimg source is at rest in 
a completely stationary homogeneous fluid of infinite extent. The 
origin of coordinates is fixed in the source, and the conventional wave 
equation applies. The velocity potential produced hy the source is 
readily expressed as the double integral of the velocity potentials for 
all reduced wave numbers, and in the s (streamwise) and y (lateral) 
directions. This corresponds to the decomposition of the source poten- 
tial into an infinite set of plane waves (and exponential disturbances). 

The source and some portion of the surrounding fluid at rest must 
next be confined within a jet. From coordinates fixed in the source 
this jet is seen as shear layers flowing past on each side, and outside 
the shear layers the ambient fluid flowing by at a fixed velocity. (See 
Fig, 1,) To accomplish this insertion of the source into the jet one 
must add to its original velocity potential the potentials for upward- 
moving and downward-moving waves reflected off the shear layers. These 
two reflected wave amplitudes are as yet unknown. 

The Shear Layers and Ambient Air 

In the shear layers the conventional wave equation does not 
apply. The correct partial differential equation is derived, and by 
assuming periodic solutions in the streamwise and laterad directions 
an ordinary differential equation is obtained in the coordinate normal 
to the shear layers. The ordinary differential equation is solved by 
power series expansions about the singular point, or about other points, 
and two Independent solutions are obtained (in each layer) with ampli- 
tudes as yet unknown. 

In the ambient air the conventional wave equation applies again 
(for coordinates fixed in the ambient fluid). Only outward moving 
waves need be considered, so (in each ambient region) one solution of 
unknown amplitude appears. 

There are now eight unknown amplitudes to be determined and four 
botmdaries between fluid layers. Across each of these boundaries pres- 
sure and displacement must be continuous, yielding the necessary eight 
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Fig. 1 COORDINATE SYSTEM FIXED IN SOURCE 
AND CORRESPONDING VELOCITY PROFILE 







equations to make the system determinate, A double integral (over 
and Kg) for velocity potential in the ambient air is then obtained as 
a function of source strength and frequency. 

The Far Field 

Only in the far field can this double integral be evaluated 
easily. There the integrand generally consists of a slowly varying func- 
tion multiplying a function which oscillates rapidly about zero. Such 
functions interact weakly, and significant contributions to the integral 
occur only when (a) the slowly varying function becomes rapidly varying 
(i.e. in the neighborhood of singular points) or (b) when the rapidly 
oscillating function ceases to be rapidly oscillating (l,e, in the neigh- 
borhood of stationary phase points) . Evaluation of the Integral yields 
the pressure in the far field for coordinates fixed in the source. 

Transient Sources and Retarded Coordinates (the far field ) 

A more useful result would be the far-field mean-square pres- 
sure produced by sources in a localized region, say Immediately behind a 
jet nozzle. To obtain this we consider a sequence of transient sources, 
each originating at the same point relative to the ambient air or nozzle. 
As one source disappears after drifting downstream with the fluid, a new 
one (with random phase relative to the first) appears at the upstream 
point. The mathematical analysis of this process (Ref. 13) is rather 
tedious. However the practical application conforms to a simple rule if 
it is assumed (as in the present analysis) that the transient sources 
have a lifetime of many cycles. The rule is that the mean-square pressure 
in source coordinates should first be formally transfoinned to retarded 
coordinates, fixed in the nozzle, (See Fig. 2.) The result must then be 
multiplied by (l + MXg/fi) , where M is the source Mach number relative 
to the ambient air, Xg is the streamwise distance and R the radius to the 
far-field observation point.* (The rule holds also for the supersonic 

*The factor 1 1 + MXg/R 1 is often written ( 1 + M cos© | where 9 is 

the angular position of the far-field observation point relative to the 
jet axis, 0=0 being measured upstream , 0 = rr downstream . 

The travel distance of each transient source and the thickness of 
the jet are of course negligible compared to R. 
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Nozzle 


Fig. 2 COORDINATE SYSTEM FIXED RELATIVE TO 
NOZZLE (RETARDED COORDINATES) 
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case if the 1 1 + Mx^/R | factors aie enclosed hy absolute magnitude 
signs.) 

Solutions 

The methods just described have been amply illustrated in 
1 13 

previous reports * , Power series solutions for sources in two- 

dimensional shear-layer jets and in circular-cylindrical shear-layer 
jets are given in Ref. 1 and are not repeated here. These previously- 
derived solutions were programmed for machine computation and are the 
basis for calculating all shear-layer jet results presented in the 
first section of Part 1 of the present report. 

The partial differential equation for the shear layer is 
(37) of NASA CR-2390* This should provide a good first approx- 
imation for cold jets of low supersonic Mach number, where temperature 
gradient effects are considerably smaller than velocity gradient 
effects. For high supersonic Mach numbers and heated jets it may be 
necessary to use Eq. (lO?) of CR-2390 which includes temperature 
gradients. 

The Near Field 

Near-field calculations are much more difficult than far- 
field calculations for two reasons. Stationary phase approximations 
cannot be used, and the transient sources cannot be assumed to pass 
through many cycles. To get a realistic picture of the near field, 
transient sources must generally be required to have veiry short 
lifetimes. This will be discussed further in Part 2. 
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Part 1 PEAKS IN THE FAR -FIELD DIREGTIVITY PATTERN FOR SUPERSONIC JETS 
Introduction and Examples 

In one of our first investigations of supersonic jets (see Ref. 14) 
the exajnple chosen was for Mach number 3i the jet was two-dimensional and 
the source was located at the center of the jet. The velocity profile 
was a linear variation from maximum velocity at the center to zero velo- 
city at the edge. In this case the mean-square pressure plotted against 
angle in the far field showed zero at the Mach angle, but showed a high 
finite peak in the "subsonic approach region" and a much lower peak in the 
"supersonic approach region". Both of these peaks were near the Mach 
angle. (The source Mach number of 3 projected onto any angular ray gives 
an "approach" Mach number of 3 or less for that ray. When the approach 
Mach number is less than unity we say the far field point is in the 
subsonic approach region. When the approach Mach number is greater than 
unity the far field point is said to be in the supersonic approach range 
of angles.) 

This result seemed reasonable. The cruder representation of the 
jet by stationary fluid with "moving sources" passing through it shows an 
infinite (and unrealistic) peak at the Mach angle. The requirement that 
jet fluid must move with the sources (as in the present analysis) ap- 
parently produces a more realistic result. 

However the consideration of an off-center source in the jet de- 
scribed above (a more difficult problem) produced an unexpected result. 
Pig. 3 shows the effect of locating the source one-quarter of the way 
from the centerline to the edge of the jet. The source, drifting with the 
local fluid in the shear layer, now travels at a Mach number of 2.25 • 

The expected and comparatively broad peaks in the subsonic and supersonic 
approach regions are much lower than before and now there appears a 
"spike" confined to less than one degree of angular range but hundreds of 
times higher than the other peaks. In an ordinary survey of the far 
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field with a five-degree angle interval, or even a one-degree interval, 
such an effect might he completely overlooked. It is not yet known 
whether the area under such a spike would contribute significantly to the 
final integrated results. 

Studies were also made of two-dimensional jets at a Mach number of 
1,5 • In Pig, 4 mean- square pressure is plotted against far-field angle 
for a jet whose velocity profile is linear, giving M = 1,5 at the center 
and M = 0 at the edge. The source is at the center, A high peak appears 
in the subsonic approach region as in the M = 3 case, and the supersonic 
approach peak is negligible (see Fig, 5)» Pon a second case we let the 
jet have uniform velocity over the central half of the jet thickness, then 
velocity decreasing linearly to zero at the edge. The source is located 
halfway between the center and the edge, and so its Mach number is 1,5 as 
in the preceding case. However the subsonic approach peak is now greatly 
reduced (see Fig, and a supersonic approach peak of similar magnitude 
appears (see Fig, 5)* Such effects are not easily explained. In a third 
example the source is located in the shear layer, three-quarters of the 
way from the centerline to the edge. This source, drifting with the 
fluid, travels at M = 0,75 • There is of course no supersonic approach 
region based on source Mach number, and the mean-square pressure plot 
somewhat resembles the subsonic jet type as might be expected (see Fig, 6) 

We now turn our attention to circular cylindrical jets with M = 1,5 • 
In each case the jet has a central core traveling at uniform velocity 
(corresponding to M = 1,5), and an annular region in which the velocity 
varies linearly from the central value to zero at the outer edge of the 
jet. The radius of the central region is Tq , and the radius at which 
the source is located is r^ , The mean-square pressure is an average 
value obtained by integrating around the jet. In Figs, ? and 8 results 
are shown for a variety of values of core radius, the source being located 
either at the center or at the edge of the core. Thus the source, drift- 
ing with the fluid always travels at M = 1.5* In Fig, 7 the subsonic 
approach region is shown. The subsonic approach peak is highest for the 
smallest central, core, and decreases rapidly as the core size is increased 




Pig. 4 EFFECT OF VELOCITY PROFILE ON FAR -FIELD MEAN-SQUARE PRESSURE 
FOR SOURCES IN TWO-DIMENSIONAL SUPERSONIC JET 
(Subsonic Approach Angles, M^ = M^. = 1,5> Sp “ 0 . 2 ) 

o J 
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Pig, 6 FAR -FIELD MEAN-SQUARE PRESSURE FROM SOURCES IN 

SHEAR LAYER OF TWO-DIMENSIONAL SUPERSONIC JET 

(M^ = 1.5, M = 0.75, = 0.2) 
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131i 8 - Mach angle for M. 
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FAR-PIELD MEAN-SQUARE PRESSURE FROM SOURCES IN THE 
CENTRAL REGION OP A SUPERSONIC CIRCULAR JET 
( Subsonic Approach Angles, M = M. = 1.5, S_ = 0,: 










In fig, 8 the supersonic approach region is shown. It is necessary at 
this point to explain that a source at radius r and angle ^ = 0 in 
the jet is constructed hy the Fourier series representation of a delta- 
fvinction. The first term is a uniform distribution of source strength 
on the circle of radius r^ , The next term is a simple cosine term, 
the next a cosine term, etc,, each with the appropriate coefficient 
to represent a delta-function. Since the wavelength created by the pul- 
sating source is generally much greater than 2r^ we expected the zero 
harmonic to behave much like a simple source, the first harmonic like a 
dipole, etc. Each higher harmonic would then have a much smaller con- 
tribution to the mean-square pressure than the preceding one. In many 
cases this turned out to be correct, and convergence was rapid. However 
Fig. 8 shows for r^ = r^ = 0,5 a very high peak due largely to the 

third harmonic*. This is a rather surprising result. Presumably it is 

produced by some sort of quasi-resonant effects associated with commun- 
ication around the jet. Such effects would not appear in the two- 
dimensional examples. However the "spikes" occurring in two-dimensional 
examples also occur for the circular jet, though perhaps more rarely. 

For r„ = r = 0,^5 > Sm =0,2 , M. = 1,5 (see Fig. 9) a spike was en- 

Us 1 J 

countered in the supersonic approach region. This very high peak was 
contained within an angular range of one degree or less. 

In the examples just considered the source was either at the center 
or at the edge of the core, traveling at a Mach number of 1.5* In Fig, 

10 we show for several core radii the effect of locating a source in the 
shear layer just far enough out so that it travels at M = 0,7. For 

jets with larger cores the sources will then be nearer the outer edge 

of the jet, but all sources travel at the same Mach number. It is inter- 
esting that the peaks vaiy widely in position and magnitude, showing the 
extreme importance of jet velocity profile in modifying the radiation 
pattern of sources convected at a specified velocity. It is also of 
interest that the mean-square pressures are nearly zero at the Mach angle 


*Where they are of interest dominant harmonics are noted on the curves 
by a 0, 1, 2, 3 etc, enclosed in a circle. 
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Fig, 9 "SPIKES” IN FAR-FIELD MEAN-SQOARE PRESSURE FROM SOURCES 
IN THE CENTRAL REGION OF A SUPERSONIC CIRCULAR JET 

(M^. = Mg = 1 . 5 , = 0.45) 
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corresponding to the jet Mach numher even though the source velocity 
is subsonic. 

Figures 11 and 12 show the radiation patterns when the sources are 
so placed in the shear layer that Mg = 0,5 and 0,3 . The mean-square 
pressures are much lower in general, because of the lower convection 
speeds. The extreme effects of velocity profile are still apparent, etnd 
some double peaks now appear. For angles less than 131.8° mean-square 
pressures are not negligible, but are not of special interest here. 

In view of the unexpected and rather perplexing results of the 
supersonic jet analysis it seemed desirable to get a check on the 
correctness of the basic formulation of the problem and on the computa- 
tions, The following two sections are devoted to the Independent study 
of simplified models of two-dimensional and circular jets. 


The Three Layer Model of the Two-Dimensional Je t 

This case is particularly valuable in giving an understanding of 
spikes because it can be studied analytically to a large extent and 
the necessary numerical calculations can be made with a small computer 
(e.g, the HP-35). 

The model is illustrated in Fig, 13 and the analysis follows the 
general procedure outlined in the introduction. The absence of a shear 
layer of course simplifies that procedure. 

The velocity potential in any region, (^, or satisfies 

the simple wave equation 

V V = ^ !!ltt (1) 

c 

for coordinates at rest in the fluid of that region. Only outgoing 
waves need be considered in the ambient air, and for simplicity the 
analysis given here treats only ordinary waves (i,e, < 1, subsonic 

approach angles), though note is made of the changes required for the 
reversed wave case > 1, supersonic approach angles). For simplicity. 


20 



Ap 

„„ „ 2.2 2 2 
80 ^ 2it a p <ti 



Fig. 11 FAR -FIELD MEAN-SQUARE PRESSURE FROM SOURCES AT M =0.5 

s 

IN SHEAR LAYER OF SUPERSONIC CIRCULAR JET 

(Mj = 1.5, = 0.2) 


21 


131*8 - Mach angle for M. 
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Fig. 12 FAB -FIELD MEAN-SQUARE PRESSURE FROM SOURCES AT M =0.3 

s 

IN SHEAR LAYER OF SUPERSONIC CIRCULAR JET 

(Mj = 1.5, = 0.2) 
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also, we consider a symmetric problem (source at z and image source at 

s 

-z ) and an anti- symmetric protlem (source at z and "sink", l,e. 180° 
s s 

out-of -phase source, at -z ) and write equations only for z > 0, The 

solution for the single soiarce at z = z is then the superposition of 

s 

a symmetric solution and an anti- symmetric solution, each of half -strength. 

For this example the source is situated in the outer layer. The 
potentials in the three regions (Fig. 13t) may then he expressed as 


= Ao»' R.P. R 


A0, = 


A0, = 


exp{io;'(K^[xQ-(MQ-M^)ct>K2y-ct+z^l+(MQ-M^)K^]^-K^)} 
t exp|io,' (K^[xQ-(MQ-Mpct>K2y-ct-zVtl+(MQ-M^)K^]^-K^)}j 
Aai* R.P. [exp|iw*[K^x+K2y-ct+ |z-zj^^]}/ 

+ B exp|iw*[^K^XtK2y-ct+2^/l-K^] 

+ D exp|lu»*[K^x:+K2y-ct-z7l-K^] } j 

A«' R.P. S exp|la>*[K^(x2+M^ct)+K2y-ct+(z-a2)7( 1-M^K^)^-K^] J- 


( 2 ) 


( 3 ) 

W 


where R.P. denotes "real part of" and the complex notation is convenient. 
In these equations A is the source strength coefficient [A = (source 

p 

strength )/8rf, the source strength being the maximum volume introduced 
per unit time]]; B, D, R, S are unknown complex amplitude coefficients 
to he determined. In the first term is the component potential due 

to the point source; the other two terms are upward-moving and downward- 
moving reflected waves. In the upper sign applies to the symmetric 

problem described above cind the lower sign applies to the anti-symmetric 
problem, thus is chosen to satisfy required conditions of symmetry 

or anti- symmetry at z = 0 (zero normal velocity, , in the 

symmetric case; zero pressure, -/oA0q^ , in the anti-symmetric case}, 
and Kg are reduced wave numbers (K^ = k^/u' , Kg = ^ 2 /"' ’ 
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w* = w/c). K = A ^2 component potentials, of 

course, the total potential in each region being the integral over all 
positive and of the component potentials shown, (in the reversed 
wave case, the term involving the square-root in the exponential of A0^ 
changes sign, so the exponent becomes 

{lo)' |K^(x2+M^ct)+K2y-ct-(z-a2)s/(l-M^K^)^-K^ j | .) 


The boundary conditions require matching at each interface the 
pressures (on the two sides of the interface) and the displacements. 
These four conditions determine the four simplltude coefficients B, D, R, 
S. The last of these is particularly of interest in describing the dis- 
turbance which is transmitted outside the jet: 


2 [] 


S = 


N 


where 


- e“^^(e^V2-l) ^ 

[Jl, = (e-V^-^l) + e-%2-1)) 


D 


and where 


f = cu'(z^-a^)7l^ 


K 


= w'(a2-a^)7l^ 


a = «u'a^7(l+[MQ-MjK^)^-K" 

_ 7(i+[Mq-mJkP^-k^ 


= 


Mo - 




(5) 


( 6 ) 


(7) 


( 8 ) 
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Again, the upper sign applies to the symmetric solution and the lower 

sign to the anti-symmetric solution, (For the reversed wave case th e 

corresp onding soluti ons are obtained by substituting -I 

for +y( wherever it occurs.) 

Application of the far-field analysis (using the principle of sta- 

1 12 l6 

tionary phase) is described in previous reports’*"’ ’ . The critical 

K- values, which alone determine S in the far field, are, in retarded 
coordinates. 


^2^ 

l+M^Xg/R 

y/H 

l+M^Xg/R 


(9) 


/ 2 2 2 

where X 2 »y,z (and R = +y +z ) are the retarded coordinates fixed 

in the nozzle (which is assumed to be stationary in the ambient air). 

From the far-field potential for a permanent source of frequency w 

(the case so far considered), the far-field potential for a transient 

source follows by Fourier integral construction. For the transient 

source sequence described in the Introduction and used in previous 
1 13 1^ 

studies'^’ the mean-square pressure in the far field, assuming each 

13 

transient source to go through many cycles of frequency , is 

— _ ^ S'^(tt,Q) + S"^(wq) 

^ R^ R^ I l+M^Xg/R I ^ 


where S = S* + iS" . For convenience the subscript 0 on the gener- 
ating frequency is dropped in the remainder of this analysis. 

Since the mean-square pressure in the far field is proportional to 
the square of the absolute value of the transmission factor S , evalu- 
ated at the critical values of K. , K_ (Eqs, (9))» attention will be 

I —12 ^ ^ 

focused primarily on 1 Sj • Results were considered only in the plane 
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y = 0 where Kg = 0 and K = Ik^I . We define 0 = 0 in the upstream 
direction so that 0 = cos"^(xg/R) . 

When Eqs. (5) and (lO) were investigated for several numerical ex- 
amples, spikes were found only when K was large (K > l) and when 

5 that is, spikes occur in this example in regions 

of the far field where the critical K-values correspond to exponential 
components of the source hut to true wave disturbances in the central, 
region of the jet. For this region, if we write f = if* , <r- icr* , 

= -i^^* , /tg = -i^' (so that f ' , cr* , t /tg* 


is|2^ 4 

where, in the symmetric case, 

N = coshf* cos a - /tg* sinhf sin a 

= sinh<r* cos a - P- 2' cosh a* sin a 

Dg = cosho* COSO - /tig* sinha* sino 

with 

a* = <l»*(ag-aJ^)^^^-l 

f* = o.*(z^-a^)J?^ 


( 11 ) 


( 12 ) 


(13) 


and, in the anti-symmetric case. 
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(l'») 

In this form it will he noticed that when K is laxge, ^ ' is small, 

Thu^ for large K, is normally the dominant term of the denominator. 
However it is oscillatory and zeros of can occur in the symmetric 
case when 


N = coshf* sina ^inh^' coso 

= sinhff* sinaf/tg* cosho-* cos a 
Dg = cosher* sina+ /ig* sinha* cos a 


tanher* = /* 2 * tana (l5a) 

and in the antl-s 3 rmmetric case when 

tano tanher’ = - /^2*' (l5h) 

It is when the dominant term of the denominator drops to zero that sharp 
peaJks occur in the far-'field mean-square pressure distribution. For very 
large K (which according to Eq, (9) occur for supersonic values of 
near the corresponding Mach cone, i,e, near 0 = cos“^(-1/m^) ) there can 
he multiple solutions of Eqs,(l5), and the pressure peaJes are very abrupt 
in 0, becoming spikes. The magnitudes of the spikes are, for S » 

S = f/sin0 (l6) 


and for the corresponding mean-square pressure 


A 2 „2 „2 

Ap R _ f 

2 2 2 2 5 

2n A p <0 1+M^cos© 


where, in the symmetric problem, 

coshf cos a - pu' sinh^* sina 

f = f £ 

® cosho-' COSO - sinho-' sina 


( 1 ?) 


( 18 ) 
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and in the anti-symmet-ric problem 


coshi* sina + /i.„* sinhf' cos<* 

f = f = ^ ; (19) 

^ cosher* slno + /tg* sinher* cos a 

Two source locations are especially of interest. If the source is 
at the outer edge of the outer layer, then i ' = a-* and = 1 • 

S cL 

? 16 

Note that Ap then corresponds to that of the modified moving source 

If the source is at the inner edge of the outer layer, then {* = 0 and 
(using Eqs, (15)) f = f = coshtr* . The mean-square pressure in this 
case is much larger than that for the modified moving source (and ap- 
proaches infinite value at the Mach cone for supersonic M^). 

Note also, however, that there axe zeros in the symmetric problem 
contribution to the fax-field pressure distribution when 

tanhf* tana = 1 /^* 2 ' (20a) 

and in the anti-symmetric contribution when 

tana = tanhf* (20b) 

For very large K these zeros are not fax removed (in 0) from the spikes 
located according to Eqs, (l5a,b). 

Further examination of Eq, (5) reveals additional information for 
other parts of the fax field where spikes were not found, 

(a) When > 1 and > (1 +[Mq-M^]K^)^ , then neither 

term of the denominator of Isl^ is oscillatory and no spikes can occur, 

(b) When < 1 and < (1 +[Mq-M^3k^)^ , both terms of the 
denominator are oscillatory but neither is dominant (i.e, I/a^I =0(1) ), 
Hence, although pressure peaks may occur, their magnitude is more limited. 
Also, since the variation of the critical K-value with 9 (location in the 
far field) is not rapid for those 0 where K < 1 , there axe no real 
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spikes • Pressiire peaks occiir singly and with moderate slope (in plots 
of Ap vs, 0), 

This example illustrates the type of problem encountered when the 
supersonic shear-layer jets were Investigated, For the shear-layer jets 
abrupt spikes sometimes occurred in the numerical solution, as illus- 
trated earlier. For all cases checked, the spikes were associated with 
the sudden vanishing of one term of the denominator, a term which dom- 
inated the denominator at nearby 0's, These spikes existed over less 
than a degree in 0, so that they could easily be overlooked in any but 
an extremely fine survey of the far-field. The question then arises 
whether it is possible to predict that section (or sections) of the fax 
field where spikes might be expected and whether one can know when all 
of the spikes have been located for a given problem. The very simple 
example being considered here, where the jet shear layers are replaced 
by two constant-velocity layers on either side of a central layer of 
higher constant velocity, shows some of the difficulties of answering 
this question. 

For this example, the three-layer model of the two-dimensional jet, 

the following generalities probably can be made: 

2 

(a) Spikes occur only when >1, (This appears to be 

true also for shear-layer jets,) Now, in general. 


K. = 


COS0 


1 1+M^cos0 


so that if is subsonic (M^ < l) 


1 1 

rr > Q >0 corresponds to - Yljif < 


(21a) 


and if > 1 


rr 


to 

M^+1 ^^1 

to 

^ > Ki > j 


(21b) 
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(See Fig. l4-a.) 


(Id) The possibility of spikes requires that 
he real. If (Mq-M^^) < 1 this happens for 


and if (Mq-M^) >1 it requires 


— oo 


< K. < - 


1 ^ “ (Mq-M^)-I ’ 


(Mq-M^)^ 


< K. < oo 

1 


(22h) 


(See Fig. l4-h.) 

Thus, if (Mq-M^) is subsonic and is subsonic, there is no 

overlapping of the K^-ranges of Eq. (21a) and Eq. (22a) such that > 1, 

and no spikes can occur. 


If M> < 1 and (M„-M.) > 1 there will be an overlapping of the 
1 ' 01 ' 2 
K^-ranges of Eq, (21a) and (22b) in a region where >1 if 

1 

1-M^ (Mq-M^)-I 

This requires 

Mq >2 

(Mq is the jet center-region Mach number.) Note that the overlapping 

Kj^-range is finite, so there can be relatively few peaks, if any. Also, 

although > 1 , it is not excessively large and the magnitude of the 

peaks is therefore limited. Since the phase change of a (Eq. (8)) is 

limited in the overlapping K^-range, it is not difficult to numerically 

locate all of the peaks in this case, using Eqs. (15)« 

If M. > 1 and (M^-M. ) < 1 there will again be overlapping of 
1 1 2 
the pertinent K^-ranges (Eq, (21b) and Eq, (22a)) when > 1 if 



VoJ 

^0 





Again this leads to the requirement Mq > 2 . Again the peaks are 
limited in number and somewhat limited in magnitude. 

Finally, if both > 1 and > 1 there is overlapping of 

the Kj^-ranges of Eq, (21b) and Eq, (22b) for 

K. < and K. > — ^ 

(Mq-M^)-I ^ M^-1 

(Of course, if (M^-M^) > 2 , we are interested only in < -1 on the 

negative side, sind if M. > 2 we are interested only in K> > 1 on the 

1 2 ^ 

positive side, in order to stay within the >1 criterion.) Since 

Ik^I is unlimited in this case, the number of possible peaks is infinite 

and the magnitude of the spikes can get very large when |K^| — oo (that 

is, when 0 -*■ cos (-1 /m^) , the Mach angle corresponding to the source 
velocity). Note also that the spacing (in 0) between the spikes dimin- 
ishes continually as | K^| °o , 

We repeat that all of our calculations have been made for the plane 
y = 0 in this two-dimensional problem, so that Kg = 0 , K = 1 K^| 

(from Eqs. (9)). If Kg is not zero, that is y 0 , the foregoing 
conclusions probably can be extended to the appropriate overlapping re- 
gions of the K^ , Kg plane (see Fig, 14). The criteria for locating 
spikes are less simply written in such cases, but appear to be straight- 
forward enough. However, if the shear layer of the jet were to be 
approximated by multiple layers, each introducing phase-changing quan- 
tities similar to a , o- used above, it will be foreseen that the criteria 
for the possible existence of spikes become very complicated indeed. 

Analysis similar to the preceding has been carried out for the 
three-layer two-dimensional jet when the source is located in the central 
high-velocity layer of the jet. In this case no spikes were found in the 
far-fleld pressure distributions of any of the examples investigated. 

This result appears to follow from the fact that, although the denominator 
of the transmission factor |SP is the sum of two oscillatory terms, 
neither one is truly dominant in any K^-range (that is, in any 0-range). 
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The Circular-Cylindrical Jet Having Two Constant-Velocity Layers 

This model consists of a cylindrical central region having constant 
high velocity, surroimded hy an annular layer having smaller constant 
velocity. Viewed in any section containing the axis of the jet, the vel- 
ocity distribution appears the same as previously sketched for the corres- 
ponding two-dimensional model. Fig, 15 shows the velocity distribution 
and the parameters used to describe the two-layered circular cylindrical, 
jet. 

This model, like the previous one, has also the advantage that an 
analytical solution can be obtained. Because of the cylindrical func- 
tions involved, however, and because an off-center source is represented 
mathematically as a series of such functions, numerical calculations for 
specific examples are considerably lengthier than for the two-dimensional 
case, A survey of the far-field pressure distribution for any example 
would be formidable without the aid of a programmed computer. 

In cylindrical coordinates the velocity potential for a source 
which is at rest with respect to the immediately surrounding fluid and 

which is located at t , f may be written 

s s 


= ttA 


oo 

I E.P. y . 

'' -oo 


n=0 


(r > r ) 
^ s^ 


or 


(23) 


OO 

= rfA«' ^ cos[^n(</^-</^^)]| R.P. ^ e^“ • 

n=0 -oo 


. J^(r«'/^) dX 5 


(r < r ) 
^ s^ 


(a is the source strength factor previously defined; m* = w/c where c 
is the speed of sound; X = k/u* where k is wave number in the x- 
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direction; «q = 1, = 2 for n = 1,2,3, ••• ; J^( ), ) axe 

cylindrical functions, basic solutions of the Bessel equation.) For 
convenience we write the source potential as 

OO 

oo ^ 

= R.P. I /0s (K) dK (24) 

® n=0 J 

-OO 

where component source potential for one wave number pair 

(K,n) • (in what follows we will write only component potentials for the 
various regions, it being understood that the total potential is obtained 
from the component by integrating over all K (positive and negative) and 
summing over n (zero plus all positive integers),) 

In each of the three regions shown in Fig, 15 the velocity potential 
satisfies the wave equation (in cylindrical coordinates) for coordinates 
at rest in the fluid of that region. With the source located in the 
outer annulus of the jet we can write 

(a) for region (O < r < r^) (See Fig. 15) : 

= ttAw' cos[n(v^-<^g)] R.P.[r^ exp|iw*[KjXQ-(MQ-M^)cti-ct]| « 

• Jn(r“’V(l+[Mo-M^]K)^-K^) ] (25) 

(b) for region (r^ < r < r^): 

^ + ttAui* • cosrn(^-«(/^) 1 R.P. exp|iw*(Kx-ct)| • 

^n Sn n s -■ l . 

* C H(2)(r«'yi-K^) ]] 

(26) 
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(c) for region (r > r^)i 

= TfA«* cos[n(^-V'g)3 K.P. & exp|i«' [K(x2+M^ct)-ct]} . 

. H^^^(r„'v(l-M^K)^-K^)^ (27) 

(When (1-M^K)^ >K^ and K > (i/M^) , (rw*y(l-M^K)^-K^) is 

replaced by H^^^(rw'vtl”M^K)^-K^) ; this is the "reversed wave" re- 
12 ^ 

gion , and the change to Hankel functions of the second kind preserves 
out-going waves as required for the solution outside the jet.) 

Satisfaction of boundary conditions of equal pressure and equal 
displacement on either side of the boundaries at r = r^ and at r = r^ 
permits evaluation of the complex amplitude coefficients R^, E^, G^, S^, 
Using the following convenient, definitions 


'0 - "’^0^ 
cr^ = 

^ = «'rQ7(l+[MQ-M^]K)^-K^ (28) 

77 = «'r^V(l-M^K)^-K^ 

Vu-m^k)^-k^ 

= 

(l-M^K)^X^ 

y(l+[MQ-MjK)2-K^ 

^2 " 

(i+[mo-mJk)Vi-k^ 
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( 29 ) 


the transmission coefficient for the cylindrical layered jet is 


where 




s = 


'2n 


^ TTcr^(l-M^K) 


[] 


■'ll 


[ J 


[ J 


N. = -/*2 


*-n 


= Jn(>>) •!„ (o-o) - J*2 J„ (X) J„(o-o) (30) 

2n 

Hd„ - kO') - "2 ' 

-[j„(x) -'‘2 '■k('o)] • 


(Note again that, in the reversed wave region \"h) i'ts deriv- 

ative (rj) are replaced by (•»?) ,) 

The far-field analysis, the change to a transient source sequence 

(instead of a permanent source) and to retarded coordinates then follow 

1 14 

in the manner indicated in previous studies ’ , As shown in Ref. 1 

(Eq, (95)) t the average mean-square pressure in the far field due to a 

ring of random-phase transient sources (of many cycles duration) located 

at r = r is 
s 


Ap^ R^ 
2.2 2 2 


Ztt A p o 1+M^cos0 


-:r y € (S *^+S "^) 
|3 ^ n' n n ' 


n=0 


(31) 


39 


where S *,S " axe evaluated at the critical value of K which for a 
n n 

point ^2’^ ^ field is 

K = cos0/(l+M^cos0) ( 32 ) 

(0 = cos’^Xg/R , R = Vx^^+r^ ; = S^' + iS^" ). 

1—12 2 2 

To discover peaks we would examine S = S ' + S " . For the 

^ n n n 

region > 1 , (1 +[Mq-M^]k)^ > Y? , for example, Eq. ( 30 ) leads to 



2 


4 





(33) 


where 


N = I (<r •) fj (x) K - Uo' j K (o-n')l 

n n^s^Ln^' n^O'^ ^2 n^'^ n^O'^J 

- K («r ') fj (X) I '(<r-*) - fi' J '(X) I (o-^')! 
n^ s L n^ n ' 0 ^^2 n ^ n^ 0 


= 

■*-n 


"'n'('-o') -'‘ 2 ' V'o’’] • 

• -'‘1’ ■^n'M InK')] 

-'* 2 ' Jn'O-) • 


(34) 


\ - [V’') - >‘2 VM V<'o’’] • 

• [V’5 \'^’i"> - '‘i’ \'M i„(»i’)J 

- -'‘ 2 ’ y°b')] • 

• [V’>) - '‘l' V'l')] 
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with = io-Q*. O-^ = = ~W 

that Oq', ®‘g'» ^2* real. (Primes on the Bessel 

functions ^n*^ ^n'^ indicate a derivative 

with respect to the argument; i.e. ~ 

L. J I 

1 

Although D. and both are oscillatory, it is not obvious that 
in ^n 

either one is dominant in some K-range (nor that it at the same time 
has one or more zeros within that K-range). Such a situation may occur, 
hut it is prohahly much less frequent than in the corresponding two- 
dimensional case studied previously. 

In the limited number of numerical examples investigated, no spikes 
in the mean-square pressure distributions were found. Even for these 
examples, however, it is not clear that spikes do not occur at all. The 
far-field survey may simply have over-looked them by using too-large 0- 
increments. (Preliminary inspection of computed results indicates that 
this is probably the case, the calculations showing denominator factors 
that oscillate very rapidly and with wide amplitude variation as © 
changes, in some instances.) These examples did show important harmonic 
effects. That is, for ©-values where K is large, the primary contri- 
bution to mean-square pressure, from Eq. (3l)i did not come from the 
n = 0 source- strength harmonic alone; contributions of approximately 
equal magnitude often were made by the first several harmonics. 


General Discussion of Peaks 

The preceding sections indicate that the peaks or "spikes" appearing 
in plots of mean-square pressure versus angle for supersonic jets are not 
the result of numerical errors or computational difficulties, nor are 
they the result of convergence problems in series solutions for the shear 
layer. Instead they constitute a real part of the theoretical results. 
(Also, the lower and broader double peaks produced by harmonic effects 
in the circular jet appear to be real.) 

It has already been suggested that the spikes axe quasi-resonant 
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effects, and in Appendix A a drastically simplified model is shown to 
exhibit both resonant and quasi-resonant effects. (A denominator passing 
through zero to give an infinite peak corresponds to resonance. The 
dominant part of a denominator passing through zero to give a finite peak 
corresponds to quasi-resonance.) 

The two-dimensional layer model, which can be studied analytically, 
offers a much better method for investigating peaks (see earlier sub- 
section). For this model only sources in the outer layers create peaks, 
and these occur only when there are e3cponential-type disturbances in 
the outer layers and true waves in the inner layer ( exponential -wave) . 

The three other cases (wave-wave, exponent ial-esqponential and wave- 
exponential.) do not produce peaks. (The disturbances outside the jet 
must of course always be true waves if they are to persist into the far 
field.) 

The peaks occur near the Mach angle for the source Mach number, 
and in some cases there may be infinitely many peaks. The observed 
frequencies are generally high, approaching Infinity at the source Mach 
angle. The peaks themselves are finite (except precisely at the source 
Mach angle in the layer model case) and always occur when the dominant 
one of two denominator terms passes through zero. 

When the sources are at the outer edges of the outer layers the 
peak values of mean-square pressure are the sajne as those produced by 

"modified moving sources" , designated "point volume displacement 

5 9 

sources by Mani and by Morfey and Tanna . 

Experimental observations in the fax field of a supersonic jet 
would probably never reveal "spikes". Such observations include the 
effects of many sources having different positions, velocities and 
frequencies, and the smoothing produced by this Inevitable "integration" 
should conceal any spikes. For the same reason a mathematical integration 
over frequency (e.g. that required to represent a transient source of 
short lifetime) should eliminate spikes in the final theoretical results. 
However this does not remove the numerical difficulties in precisely 
defining the shapes of spikes contributing to the final smoothed result. 
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Only if the areas under the spikes were negligible (compared to the area 
under the remainder of the mean-square pressure vs, 9 curve) would the 
difficulties vanish, and this has not been investigated. 

It has been noted that the spikes seem to occur at high observed 
frequencies. For example, in Figs, 3 an<i 9 the observed frequencies axe 
8,5» 6,7, and 12,7 times the generating frequency*. The simplest way to 
avoid numerical difficulties created by spikes is to exclude the higher 
ranges of observed frequency. This is not inconsistent with experimental 
approaches since the radiated energy in the far field is often recorded 
for several observed frequency bands. 

To pursue this idea further one might attempt to plot far-field 
mean-square pressure versus angle for prescribed observed frequencies 
(or observed frequency bands). However this requires the assumption of 
source intensity as a function of generating frequency. (A fixed observed 
frequency is produced by one generating frequency at one angle and by a 
different generating frequency at another angle.) Such an assumption 
might require extensive investigation. 

Heating the jet can have the effect of reducing the probability of 
encountering spikes in the pressure distribution since for fixed velocity 
the local Mach number (based on local temperature) is then reduced. This 
is illustrated in Appendix B for the simplified three-layer two-dimensional 
jet model when each layer has a constant temperature higher than ambient. 


* The Strouhal number is 0.2 based on generating frequency, local jet 
velocity (at jet center) and local jet diameter. The Strouhal numbers 
based on observed frequencies axe then 1.7, 1.3^ and 2.5^. 


CONCLUSIONS - Part 1 


Two une 3 q>ected features appear in plots of mean-square 
pressure vs. angle in the far" field for supersonic jets when 
sources axe moved off center. First, "spikes" (very high, very 
narrow peaks) are found near the Mach angle for the source for 
two-dimensional jets and, to a lesser extent, for circulax jets. 
Second, for circular jets, lower and "broader multiple peaks 
appear. These are associated with harmonic effects (communi- 
cation around the jet) . 

Neither of these features is the result of numerical errors 
or computing difficulties, nor are they the result of convergence 
pro"blems in series solutions for the shear layer. Instead they 
appeax to constitute a real paxt of the theoretical results for 
cold, low supersonic Mach number jets. However further investi- 
gation is certainly desirable to determine if the areas under the 
spikes axe large enough to provide significant contributions to 
integrated results, and to consider the effect of temperature and 
temperature gradients (see Eq. (10?) of NASA CR-2390 ) for high 
supersonic Mach numbers and heated jets. 

Examination of a simplified (two-dimensional, three-layer) 
model Indicates that high temperatures in the jet can materially 
reduce the probability of encountering spikes. This is associated 
with the decrease in local Mach number for given jet velocity and 
increased speed of sound. 

Investigation of simplified models also suggests that the 
imexpected features axe quasi-xesonant effects (where the dominant 
part of a denominator passes through zero) , and that spikes appear 
only for non-uniform velocity profiles. 
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Integration over source position, velocity and frequency 
(which occurs automatically in experimental ohseivatlons) should 
smooth these peaks. However the "spikes", hecause they may appear 
and disappear within a one-degree (or less) angular range in the 
far field, present computational difficulties. 

6. The"spikes" are apparently associated with very high ob- 

served frequencies, and at present it may he necessary to ignore 
this high frequency range in order to avoid computing problems in 
some examples. 

?• Although some of the results of Part 1 are both unexpected 

and imperfectly understood they cannot safely be ignored. Further 
study is necessary to show their true significance. 
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Part 2 NEAR-FIELD DISTURBAITCES FROM TRMSIEHT SOURCES IN SUBSONIC JETS 


Introduction 

Two factors make the near-field investigation much more difficult 
than the study of the far field. First, it is not possible to choose a 
critical (or dominant) K for each point in the near field. The at- 
tempt to do so leads to non -localized and hence unrealistic near fields. 
The disturbance produced by any one source should tend to decrease as 
the observation point is moved away from the source (upstream or down- 
stream), and this decrease does not appear for a single K, The second 
factor is that transient sources which pass through many cycles (a simple 
case which avoids a frequency integration) cannot be used because they 
also lead to a non-localized near field. Instead the sources must be 
restricted to lifetimes of the or^er of one cycle and an integration 
over frequency is required, A new assumption must be made here, the 
length of time the source exists. 

Another complication associated with the near field is that regions 
of exponential decay cannot be immediately discarded as in the far field. 
One of these regions (Region F , with true waves in the jet and expo- 
nential decay outside; see Fig, 18) may contain resonant points (singu- 
larities) which offer some difficulties in numerical calculations. 

In the following analyses we consider several cases in order of 
their difficulty. First, the stationary pulsating source in a homogeneous 
fluid of infinite extent is studied. Next, we examine the harmonic source 
in a uniform jet and, finally, the sequence of transient sources in a 
uniform jet. 

In general we look at mean-square radial displacement and mean 
energy flow. However the energy flow for the harmonic source in the jet 
is omitted because this involves a very tedious hand calculation, (Pro- 
grajnming for machine computation seemed inadvisable at this time,) For 
the transient source space-time correlations are considered also. Numer- 
ical results for the transient source case are limited, primarily because 
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e3q)erlenced programming assistance was not available to this project 


The Near Field of a Stati onary Source in a Homogeneous Fluid 
of Infinite Extent 

This solution is almost trivial, involving the well-known spherically 
symmetrical pulsating source. However we choose to discuss it in the cyl- 
indrical coordinates appropriate for a jet. 

In spherical coordinates 

0 = ^ sinp(R-ct)] (35) 

where ^ is velocity potential i A is an arbitrary amplitude, R is 
spherical radius, w = frequency, w* = cu/c, c = speed of sound and 
t - time. 


R = ^x^+y^+z^ (36) 

where x,y,z are Cartesian coordinates, or 

R = Jx^+x^ ( 37 ) 

where r is the cylindrical radius and x the axial direction. 

In cylindrical coordinates 

^ — slnru>'(.^/?+? - ct)] (38) 

To study instantaneous local energy flow, E , across the cylinder 
walls we need Ap , where Ap is the perturbation pressure, 

Ap = - Pq t with pQ - mass density of the fluid. 
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( 39 ) 


To study the radial displacement, 17 , we use 


V = f 0^ 


(40) 


The time mean of the energy flow, E, is 

o 2,2 ,2 

_ ZfiQir A o»' or 

E — 3/2 

/ 2 , 

(x +r ) 


(M) 


The time mean of the radial displacement squared, 17 , is 

..2 


— 2 2 2 
_2 2it A r*^ 

rf = 


w 


w 


2 3 

/ 2^ z-r , 2, 2^^ 

(x +r j (x +r ) 


Let X = wavelength == 2Tf/ui' , w' = w/c , Then 

.2 


^2TcY _ (Vr)-^ 

\^k] 2 2 

8(1 + 

r l_ 


2rr 

x/r , 


2 

r 


(42) 


(^3) 


In Fig, (16) we normalize to unity at x/r = 0 , and show how the 
mean-square radial displacement, Eq. (43), falls off for points "up- 
stream" or "downstream" on the fictitious cylinder walls. This is done 
for X /r -* «> and for X/r 0 . Also we show from Eq. (4l) how the 
normalized energy flow varies with x/r , The curves are symmetrical 
about x/r = 0 . 

These curves show source effects to be fairly well localized in 
this simple case, and can be used as a basis for comparison in later 
cases. 
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nor: 


Meaji Energy Flow 




Fig. l6 SOURCE EFFECTS AT FICTITIOUS CYLINDER WALL FOR A 
STATIONARY SOURCE IN A FLUID OF INFINITE EXTENT 




The Near Field of a Permanent Pulsating Source at the 
Center of a Uniform Girculax Jet 

Introduction 

Here we consider a permanent pulsating source (a harmonic 
source), on the centerline of a circular cylindrical jet. The jet vel- 
ocity is uniform throughout, and the source drifts with the fluid per- 
petually, The origin of coordinates is fixed in the source, and we 
examine disturbances on the jet boundary upstream and downstream from 
the source. Mean-square radial displacements are compared with those 
for the previously studied stationary source in a homogeneous fluid of 
infinite extent (a zero-velocity jet case). 

This problem requires fairly extensive calculations, but is still 
much simpler than the transient source problem treated later. Results 
are given for one example, obtained by hand-calculation using an HP-35 
computer. 

Development 

The origin of coordinates is fixed in the source (see Pig, 17) 
with X the streamwlse distance from the source and r the radial dis- 
tance, In this coordinate system all the fluid within the jet is at 
rest and the acoustic equation is 

\ (44) 

C 

where 0 is the velocity potential, V is the Laplacian, c is the 
speed of sound and t is time. 

The potential for a source at the origin is 

-foo 

jZlg = trA«* R.P. f H^^^(ro,'X^) (45) 

where A = (source strength)/8Tr^ , source strength being defined as 
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maximum volume per unit time, oi is the source frequency, u>* = w/c , 

k is wave number in the x-direction, K = k/w' , Hq "'(o) is the Hankel 
function of the first kind of order zero and argument a , and R.P. 
denotes "real part of". 

Allowing for reflected waves the component potential ^^(K) inside 
the jet is 

(if6) 


0(k) = ttAw* R.P. 


Iw' (Kx-ct) 





where R = R'+iR" , a complex reflection coefficient and Jq is a zero 
order Bessel function. 

Outside the jet the simple wave equation (Eq. (^^)) must be satis- 
fied for coordinates fixed in the ambient air, x^ and r . The com- 
ponent potential, > allowing for outgoing waves only, is then 


^^(K) = ttAw’ R.P. 


ia>'[K(xp+Mct)-ct] /.N y s — p 

S e 2 Hy^(ru.'V(l-KM) -K^) 


(^ 7 ) 


where S = S'+iS" , and S' and S" are in-phase and out-of-phase 
transmission coefficients. 

The total potential outside the jet is 



oo 

/ 02 ^^^ ^ 

— OO 


To determine we must find S' and S" . This is done by 

matching at the jet boundary the inside and outside values of pressure 
and of radial displacement. The pressure is Ap = ~ Pq 

Ap ~ ~ Pq ^2^ outside. The radial displacement is S inside and 

J* ^2 outside, x being held constant in the former case and x^ 

constant in the latter. We omit the details and record the result of 
applying the boundary conditions. 
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2A f i(KM-l) (r <o*y( 1-KM)^-K^) dK 

^ ^1 -1 |7km-i)2 Jl^ ] Jq’( ) - J(i-mf - K ^ ] Jq( )] 

(48) 

where 

[ ] = [r^«o*7(l-KM)^-K^J 

( ) = (r^w'Vf^) 


If we let V ~ S ^2 with Xg constant, then 


v-v/ 

2A f 

= R.P.f 

r.c J 


7(KM-1)^-K^ giw'(Kx-ct) H^l) '(rw*y(KM-l)^-K^) dX 

[(KM-1)2 ] Jq’( ) - 7(i-km)2-k^' C J Jq( )] 

(^9) 


The calculation of 7f from Eq., (49) is time consuming. First, the 
real part of the integrand must he found for four different K regions 
(see Pig. 18). Fairly extensive numerical computations must he made to 
obtain the integrand in each of the four regions, and then the numerical 
integration must he performed. For each value of x/r^ investigated the 
integrand must he altered and a new integration performed. Since sin(wt) 
and cos(««»t) terms can he removed from under the integral sign it is con- 


2 2 

venient to caiculate -r) , the time mean of 17 . 

We do not Include these tedious steps, hut simply plot the final 
results for mean-square radial displacement at the edge of the jet versus 
streamwlse position (see Fig. 19) • In the example chosen the Mach number 
(M) is 0.5f the Strouhal number (s^) is 0.2 and the ratio of wavelength 
to jet radius (x/r^) is 20 . 


Discussion 

From Fig. 19 it is clear that the mean-square radial displace- 
ments at the jet edge fall off fairly rapidly with distance from the 
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Fig. 19 


COMPARISON OF MEAN- 
JET EDGE FOR JET 
(HARMONIC SOURCES 






\j^ 

hi 

source (upstreajn or downstream) . However the rate of decrease is consid-- 
erahly less for M = 0.5 thain for the stationary source (or zero velocity 
jet). One might speculate that reflections from the edge of the jet tend 
to prevent the source effects from leaving the jet cylinder. 

It would be interesting to study the energy flow here, but this is 
a more difficult numerical problem and time does not permit. The energy 
flow should show greater asymmetzy between the upstream and downstream 
directions. (An almost imperceptible asymmetry exists in Fig, I 9 for the 
harmonic source in the jet.) 

It would ialso be interesting to study a shear layer case to compare 
with the uniform velocity jet. 

It must be remembered that the source in the jet is a continuously 
operating source (a harmonic source), not a transient source. The tran- 
sient source (or sequence of transient sources) will be considered next. 
The coordinate system can then be fixed in the jet nozzle instead of 
traveling with the source. 


The Near Field of a Sequence of Transient Sources 
at the Center of a Uniform Circular Jet 

Meaji-Square Radial Displacement 

For a permanent (harmonic) source of generating frequency cj , 
let the radial displacement of a fluid particle in the ambient air be 
ri(cj,t) • Then for a transient source of frequency Wq existing from 
t = -2 titi/wq to t = +2Tfn/a>Q (2n cycles): 





■i7(w,t) sin(2Tfnw/«Q) du» 

7^ 27 

(« ) 


(50) 


This corresponds to the Fourier integral construction of a transient 


I - 
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13 

process from all frequencies of a permanent one . 

p 

77 can then he written as a double Integral using ft as a 

dummy variable for one integration: 


V ( 


00 00 

r f it "’ 
V*)-— j j — 


t) sin(2Tfnw/«-) 77(ft,t) sin(2Trni^/«„) dft du» 


■^0 0 




(51) 


■»7(w,t) and i7(ft,t) are given by Eq, (49) of the preceding section, 
(ft and K' are substituted for w and K in the second case.) For 
brevity we write them as 


17(01, t) 


00 

R.P. f F(K,w,t) dK 

-00 

00 

R.P. S F(K',ft,t) dK' 

-00 


(52) 


Substituting Eq. (52) into Eq. (5l) gives a quadruple integral for 

P ~2 

7j (ci»Q,t) . To get the mean- square value of rj , rj , for a sequence of 

transient sources we must integrate over a sufficiently large time (t) 
and divide by 4iTn/o»Q , the time interval between the arrival (at a 

fixed point in the ambient air) of the initial impulse from one tran- 
sient source and the initial impulse from the next transient source. 




“0 

Ti^n 


00 00 

f sin(2Tmw/wp,) f sin(2Trn^^^n) 

J ~ 7 T~T^ 


00 

/ 


(“ -“ 0 ) 0 




// 

-00 -00 


J Fp(K,«,t) Fp(K',ft,t) dt dK' dK dft d« (53) 


T 

2 
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where the integration in t is performed with Xg constant. (Note 
that Xg = X - Mot ,) 

It is now convenient to make the time dependence of F explicit 
so that the time integration can he performed, 

(5^) 


F(K,w,t) = f(K,«) 


f(K,w) is defined through reference to Eq. (^9)« F(K',i2,t) and f(K',i2) 

are obtained from F(K,«,t) and f(K,w) simply hy substituting K' for 
K and for w , 


f(K,u,) = 


2A '(ru)'y(l-KM)^-K^) ® ^ 


[ ] = [r^w’7(l-KM)^-K^] 


( ) = (r^u>Vl-K‘") 

(Note that the argument of the Hankel function in the numerator contains 

r instead of r. ,) Now, denoting F = F_ + iF_ , etc,, we have 
X K X 

Fp(K,w,t) = fp(K,o») cos[«t(KM-l)] - fj(K,w) sin[wt(KM-l)] 


and so 


o OO oo 

W r sin(2THiw/o» ) r sin(2Tfn^Wf>) 


L . "O C sin(2r,na*/«»Q) r 
, 2_ 2. } 


(.2_o.^2) J 


// 


r ^ 

• fj^(K,«u) fj^(K’,il) J cos^t(KM-l)] cos|flt(K'M-l)] dt 

- -T/2 

T/2 

+ fj(K,w) fj(K',il) J sin|a>t(KM-l)] sin jftt(K'M-l)] dt dK’ dK dfl dw 

/o — I / cn \ 
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Here the integrals involving sin( ) cos( ) products have "been omitted 
since the integrands are anti-symmetrical over the range -t/ 2 to 
+T/2 and contribute nothing. 

The time integration is readily performed to give terms of the 

type 

sinj[L(KM-l) - 


[w(KM-l) - 



_w(KM- 1) + il(K'M-l)] 


¥e next consider the integration in K' and for convenience let 
T become large. The integrands aipe then slowly varying functions of 
K' multiplied by functions of the above type which oscillate rapidly 
about zero as K' varies. Such functions interact weakly, and sub- 
stantial contributions to the integral occur (a) when the rapidly 
varying function ceases to be rapidly varying (e.g. at stationary 
phase points), or (b) when the slowly varying function becomes rapidly 
varying (i.e. at singular points). This case is of the latter type, 
and so we investigate the singular points defined by 

w(KM-l) + il(K’M-l) = 0 
or 

= (wKM - « + I2)^M 

and ( 59 ) 

= (- wKM + « + i2)^M 

In the neighborhood of a singular point let K' = + S . Then 

Integrating over this neighborhood produces integrals of the type 
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= / 


sin(iiM^8) d.8 


( 60 ) 


i2MS 

which cam he evaluated as rr/iiM, rememhering that T -► °o. We then get 


o OO oo CO 

r s±n(2-[ia.u/ U q) / sin(2TmiV«o^ 

^*“0^ J 7~^ J Z2 7^ 

TT nM ^ ("-"0)0 -00 


OO 

/ 


• a 




1 x> /tt- ^ ^ / wKM-to+ii \ 

+ fj(K,«) fj( jjjp- ,fi) 


, ^ (-.r \ / -wKM+w+H 

+ fp(K.«) fp( — , Jl) 


/t, \ ^ / -wKM+c*»+i2 

~0M ’ 


dK dil do) 


(61) 


2 

In Figs. 20 and 21 the preliminaxy results of calculations for *77 

2 

from Eq. (6l) are presented. Fig, 20 shows a decrease In 77 with in- 
creasing distance radially outward from the jet. This would of course 
he expected. Fig, 21 shows that, for r/r^ = 1 amd r/r^ = 3 most of 

2 

the mean-square radial displacement (77 ) is confined to the streaimwise 

stations through which the transient source passes (a distance of ten 

jet diauneters in this case) . 

2 

For r/r^ = 1 amd r/r^ = 3 there are large fluctuations in 77 , 
It is not definitely known what causes this. These ohservation positions 
aixe very close to the source path in terms of the path length (20 ^j^)» 

~2 

One might speculate that most of the contrihutlon to 77 occurs when 
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Fig. 20 VARIATION OF MEAN-SQUARE RADIAL DISPLACEMENT WITH RADIUS 





WITH DISTANCE FROM SOURCE PATH 


the source passes "beneath the observation point. Since the phase of 
the source has not yet been made random with respect to source position 
on the p«th, some fluctuations in source effect with observation point 
might possibly be attributed to this. These fluctuations seem to disap- 
pear (or be greatly reduced) for r/r^ = 10, 

Fig, 21 also shows that the downstream disturbances exceed the up- 
stream disturbances for r/r^ = 10, Such a trend is very pronounced in 
the far field, and it is quite interesting that it shows so clearly in 
the near field (assuming these preliminary calculations to be sufficiently 
accurate) , 

The rather large streamwise extent of substantial disturbances (ten 
diameters) indicates the necessity for using transient sources of very 
short lifetime to obtain realistic results in the near field. The source 
of this example has a lifetime of two cycles. It might be desirable to 
reduce this to one cycle or even one-half cycle (though some modifica- 
tion of the equations is necessary). 

Note that a similar problem arises in the more realistic expanding 
jet. If such a model is to provide the correct far-field pressures the 
sources should move with the jet fluid. Therefore the sources should 
be transient and have very short effective lifetimes. 

For either the cylindrical or the expanding jet it may be possible 
that some simulation of convective effects can be achieved using sources 
stationary with respect to the ambient air. Such sources are in motion 

with respect to the immediately surrovinding fluid (i,e, the jet) and 

15 

so are of many types , Each type has a different relation between 
downstream and upstream radiation. It is possible that some type might 
give usable pseudo-convective effects, but this is only a conjecture at 
present, and thorough investigation would be required. 

An alternative form of Eq, (6l) which is possibly more convenient 
is obtained by introducing the observed frequency, w*, for an observer 
at rest in the ambient air, (a» is the generating frequency,) 

V = (62) 
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where M is the source Mach number. Here the subscript s may be 
s 

dropped since the jet has a uniform Mach number and the source travels 
with it. (Note that for K = = cos0/(l+MgCOs0) , «* = w/(l+MgCos0) , 

which checks the expression for w* previously used in the far field, ) 
Prom Eq., ( 62 ) 


K = («-tt)*)/Mui and dK = - dw*/Mw 

Eliminating K from Eq., ( 61 ) gives 

00 00 < 

3 


2 / 
1 ( 


U}^ f sin(2Tfnw/wQ) f sln(2Tmfya)^) C 






Mw * ^R^ MXi ’ ^ 


f,( 
- fx( 


w-w* 

Mu) 

. ) 


i2+o;* 

Mn 

. ii) 



OJ-W* 

Mw 



Mil 

, ii) 



-Jf 

CO- cu 
Mu> 

1 ") 


il+ 

Mil 

f 

dw* dSl du) 

(64) 


A triple numerical integration is still required to obtain the mean- 
square radial displacement at any point outside the jet. However it is 
possible that the response in a narrow band of observed frequencies in 
the neighborhood of w* might be obtained by a double integration. 

(Note that the for - w* must be added to that for + w* in 

such cases.) 
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Mean Energy Flow 

The rate of energy flow (E) across a unit area of a circular 
cylinder is the product of the perturbation pressure times the radial 
velocity, or 


E = 


■ '’0 


( 65 ) 


— 2 
It is E and E which we want to find now instead of 17 and 

V • (v~ f ^2 with Xg constant.) 

Since the analysis proceeds almost exactly as before we merely 
note that instead of f(K,w) and f(K',ii) we now have g^(K,«), g2(K',ii), 
where 


iw'Kx- 


gl(K,u») = 


2A«' i(KM-l) 7(1-KM)^-K^ ' (rw'7(l-KM)^-K^) 

'1 [(KM-1)^ Xi? [ ] Jq’( ) - Jq( )] 

( 66 ) 


ggCK’.il) = 


-2Ali (K'M-1)^ H^^^(rii'7(l-K'M)^-K'^) e ^ 

r^[(K'M-l)^ ^0^^ n Jq'^ ^ - y(l-K'M)^-K'^ ] J^( )] 

(67) 


(Note that ii' = 0,/c just as w' = w/c , [ ] and ( ) in Eq, (66) 

are defined in Eq, (56). The same definitions apply in Eq, (67) if K' 
is substituted for K and ii* for w' .) 
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E(«q) = 


o OO oo 

PQ^O f sin(2-nticii/wQ) f sin(2TfniV^ci>Q) 


2 „2 
TT nM 


/ sin(2Tflic<i/wQ) r 

7*2 2^ J 

0 -"o ^ 0 




OO 

/■ 


0 ^ 


[ 


/ u-y la") s ( SI) 

Mm ’ ^ ®2 r^ Mil ’ ^ 


+ a- ( MZjt w’) g. ( 11') 

®lpj'' Mm ’ ^2 r'' Mil » >* 


2r' 


+ g ( u) s: ( St) 

®ll'' Moi ’ ®2i'' MU ’ ^ 


_ g ( a>zjl w) g ( (ktsit 
Mo) ’ ^ ^2i'' MU 


, «)] 


dm* dH dm 


E(mQ) is the time mean of the energy flow across the cylinder wall per 
unit area and per unit time, oiq is the frequency for each transient 
source, which exists for just 2n cycles. Formally n should he an 
Integer, but actually correct answers are obtained for one cycle when 
n = 1/2 . w* is the frequency recorded by an observer at rest in the 
ambient air. 


Space-Time Correlations in the Near Field 
Introduction 

One part of this work was directed towards the determination 
of space-time correlations of the pressure gradient in the near field 
for con 5 )arison with in-house research at Langley Research Center. 

The mathematical formulation of this is not basically more diffi- 
cult than the treatment of energy flow just considered. However two 
problems arise. First, the introduction of pressure gradients (higher 


derivatives) gives rise to a convergence problem. This can he handled 
without much difficulty. Second, the introduction of two new parameters, 
a space interval and a time interval, tends to greatly increase the vol- 
ume of data to he obtained. 

Again we consider first the simple case of a pulsating source at 
rest in a homogeneous fluid. Then transient sources at rest in a 
homogeneous fluid are investigated. Finally, the sequence of transient 
sources at the center of a uniform-velocity circular jet is studied. 


The Stationary Source in a Homogeneous Fluid 

The velocity potential for this source in cylindrical coor- 
dinates (see Eq, (38)) is 


2ttA 

0 = sin 




2^ 2 
+r 


|w'(^^+r^ - ct)J 


(69) 


The pressure gradient in the radial (r) direction is Ap^ = - 

and if we retain only the portion associated with true sound Ap^ is 


2TfA/9f,ww’r p ^ -T 

Ap (x,t) = — - — sin|_w'(Vx'^+r'^ - ct)J 


2 , 2 

X +r 


(70) 


If we let f he a space Interval in the axial (x) direction, and 
T a time interval, then 


2TrA/)^ww'r 


sin 


[w'(y(jd-f)^+r^ - c(t+r))] 


(71) 


The correlation function is 


corr, fn, = 


+T/2 

-T/2 


(72) 
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and using the two preceding equations we get 


corr, fn, = 




COS 




2 . 2 , / 

-Hr - w '^3 




— O)-^ 


(73) 


This correlation function does not die out with increasing time 
interval (t), and so is unrealistic for the turbulent jet problem. We 
next examine the transient source which is more realistic. 

If we assume that the transient source behavior is represented by 
a portion of a sine wave (see Fig, 22ai where n is set equal to l) the 
radial pressure gradient Ap^ is (using Fourier integral representation 
as in Eq, (50)) 


(,,t) = 

^ / 2, 2\ 

-rr (x +r )c 


J“ sin 


[w' - Ct)j sln(2Trn«/wQ) dw 

7 ^ Ts 

(u) -Wq ; 




At large w the Integrand contains only periodic functions (note 

2 

to in numerator and denominator) and appears to present convergence 
problems. To avoid this we choose a smoother function (see Fig, 22h) . 

2 

Since this function (cos'^(cjQt) between -tf/2wQ and +tr/2wQ , and zero 
elsewhere) is symmetrical about t = 0 instead of anti -symmetrical, a 
cosine instead of sine is needed in the basic source potential. In 
cylindrical coordinates this is 


2ttA 

8<-- 


cos 




, 2^ 2 
X +r 



(75) 


The Fourier integral representation of the radial pressure gradient now 
becomes (for the true sound) 

^«Q^(-2TfA/>Q)r ^ (u cos[w'(-^^-(-r^ - ct)J sin(Tr‘*/2“'Q) dw 


APj,(x,t) = 


Tr(x^-^r^)c 




(76) 
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SOURCE TIME HISTORIES 


!"or t 
for 


2 < (2Tf/«o)^J F(t) = 
< (tt/2«o)^5 = 


0 elsewhere 
0 elsewhere 



2 

Smoothing the ends of the source time history produces w/« in the 
integrand for large <•> and improves the convergence. 

For a sequence of transient sources the correlation function is 

T 

2 

corr. fn. = “:J7 S Ap^C^^tCft+r) Ap^(x,t) dt (77) 

T 

" 2 

where tt/wq is the time interval between the arrival at a fixed point 
of the initial impulse from one transient source and the arrival of the 
Initial impulse from the next transient source. 

T 

[' O O O CO oo 2 

lim ^ ^ f f r sin(Tfu/2u)Q) sinCtrVSwQ) 

c^(x^+r^)[(xfC)^+r^] 0 0 T <^) 

" 2 

• cos [oj'(-^?+r^ - ct^ CO s]^* ( ^( xf f ) ^+r ^ - cit+rj ^ dt dH dw 

(78) 

The time integration can he made immediately. As T oo the SI 
integration involves a function oscillating rapidly about zero, and the 
integral is evaluated from contributions near the singular point SI = a) , 

The result is 


corr, fn. 


l6wQ^(4A^/)Q^r^) r (w/wq)^ sin^(Tr<«»/2«Q) cos(u»y/a»Q) d(w/a>Q) 

c^(x^+r^)[(x+e)^+r^] ^ [(“/"q)^-^P 


where 


y = — |~y(xhf )^+r^ - 
c ^ 



(79) 

(80) 


The integral can be evaluated (see Gradshteyn and Ryshiki Ref, 17 > 
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#3.728, 7»> p. ^10, and let their c = h+ € , €-^0). For 0 < y < rr 
we get 

I, ,2 2 4 2 

- 4rrA r r -i 

oorr. fn. - L=i”(2?-) + 2(r-Tf) oos(2y)J (81) 

For y > Tf the correlation function is identically zero, and the 
correlation function is ssmunetrical about y = 0 . Fig, 23 (from Eq., 

(8l)) shows the effect of y on the correlation function, 

A transient source at rest in a homogeneous fluid produces dis- 
turbances which are confined completely to an expanding spherical shell. 

The thickness of the shell is c times the lifetime of the source. If 
the time interval, t , is zero the space interval, f , can be chosen so 
that two points separated by a distance f in the x-direction are never 
in the expanding shell at the same time. The correlation is then iden- 
tically zero. Similar situations can arise of course when t is not 
zero. 

It must be remembered that in a sequence of transient sources the 
phase of any one source is assumed random relative to the other sources. 

The correlation function is then not affected by interactions between 
pairs of sources. 


Pressure Gradient Correlation Functions for a Sequence of Transient 
Sources at the Center of a Uniform Circular Jet 

In order to avoid convergence problems, we adopt a smoothed time 
history for the transient sources (see Fig, 22b), Since this is sym- 
metrical (instead of anti -symmetrical) about t = 0 it is necessary 
to modify the source potential 0^ by taking the imaginary rather than 
the real part of the complex potential (see Eq, (48)). 0^ then 
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Fig. 23 PRESSURE GRADIENT CORRELATION FUNCTION vs y FOR TRANSIENT 
SOURCES AT REST IN A FLUID OF INFINITE EXTENT 

(y is a function of f and t, the space and time intervals. 
See Eqs. (80) and (8l).) 
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( 82 ) 


2A r 

0, - - i.p./ 


i (KM-1) e 


iw*Kx, Iwt(KM-l) 


dK 


(KM-1)2 ] Jq'( ) - J Jg( ) 


[] = 

( ) = (r^w'^/^- 


K*- 

2~;;2 


(83) 


^ = 7(KM-1)^-K' 

The radial pressure gradient is - <1^^ 


"2Acu*oj 


i>n (".t) = — 

^rt r. 


cx 

:.P./ 


(KM-1)^ ^ H^^^'(ra,'V0) e 


iu»*Kx„ iwt(KM-l) 


dK 


(KM-1)2 [ ] Jq'( ) - ] Jq( ^ 

(8^) 


For the transient source we need ("Qf"^)* The Fourier integral 

“^rt 

representation of this is 


4«» 2 sin(Tfa»/2«^) dw 

iii (« t)--2-/ -s — —r-K 

rt Tf w(4wq -tj ) 


(85) 


To simplify this expression let 

CO 

12(2 ) = I.P. f G(K,oi,t,x ) dK 

rt ^ -00 

00 

0o (^l,t+r,x„+0 = I.P. S G(K’,ii,t+T,x,+ 0 dK' 

^rt ^ -00 

The space-time correlation function for the two points (x^jr),. (xg+f,r) 


(86) 


IS 
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(87) 


T 

2 


corr. fn. = t- ® — y'^Pj^C'^ift+TjX^+O Ap^^Cw.t.x^) dt 


TT T 
2 


corr. fn, = 


2 2 
16/>0 “o 


/ sin(Tr«*»/2wQ) T sin(Tfii/2wQ) r r 

a?") •' — -OO -00 


OO OO 


TT^ 


Cl?) •' ii(4cil^^— fl?) 


y Gj(K',il,t+T,x2+0 GjCK.w.t.Xg) dt dK' cJK dii dw 


( 88 ) 


T 

2 


In order to perform the time integration the term containing t 
is factored out, giving 


G(K,w,t,x^) = s(K,i»,x^) 


(89) 


G(K',il,t+T,X2+0 = g(K*,fl,X2) e 


mt(K'M-l) iii|^(K'M-l)+K'f/c] 


The imaginary parts of G required in the integrand are then 

Gj(K,w,t,X2) = gj(K,ci»,X2) cos[ut(KM-l)] + gp(K,<u,X2) sin[^t(KM-l)] 

(90) 

Gj(K',ii,t+T,X2+e) = g^{K\Sl,x^) cos[^t(K'M-l)] cos[i 2 (TK’M-T+K'e/c)] 

- gj(K',il,X2) sin|nt(K’M-l)] sin[fl(TK’M-T+K'f/c)] 

+ gjj(K*,X2,X2) sin(ftt(K*M-l)J cos |il(TK’M-T+K’e/c)] 

+ gp(K*,H, Xg) cos[ht(K'M-l)J sin[n(TK'M-T+K*f/c )3 


73 



In the integration range from -t/ 2 to +T/2 the anti- symmetric 
parts of the integrand can he omittedo Then 

T 
2 

Gj(K',ft,t+T,X2+0 Gj(K,a»,t,X2) dt = 

2 

gj(K,w,X2) j gj(K',il,X2) cos|ii(TK'M-T+K'f/c)] 

+ sin[S2(TK’N-T+K’i/c)J 

T 
2 

cos [wt(KM-l)] cos[ht(K'M-l2 dt 
T 

“ 2 

+ ^SpCK.w.Xg) j gp(K' Ax^) cos(ii(TK'M-r4-K'€/c)] 

- gj(K',ii,X2) sin[ii(TK'M-T+K'f/c)] 

T 
2 

• / sinE»t(KM-l)] sin[ftt(K'M-l)] dt (91) 

i 

" 2 






The time integration is simple and for T -► oo the integration in 
K’ is obtained (as previously discussed) from the neighborhood of the 
singular points • A triple integral (in ) remains to be deter- 

mined in the expression for the correlation function. As before, we 
choose to rewrite this in terms of observed frequency, <u*, noting that 
K = (w-w*)/Mw . 
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dcu* dSl dw 


(92) 


- 2kw*w (KM-1)^ 

r^[(KM-l)2 VG? J Jq'( ) - V© ] Jq( )] 

(93) 

(ft-u^)(C/Mc) - T 0 »* 


yg = (ft+W*)(f/Mc) + T« 0 * 


^ = y{^^-l){a>*^/I) + 2 («v«) - 1 
= |7(m2-1) + 2(«V‘-) - («*%^) 

KM-1 = - «*/" 
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CONCLUSIONS - Part 2 


1, A realistic pictxire of the near field must show disturbances 
dying out as the distance from a small region occupied by transient 
sources increases. The lifetime of each transient source must then 
be very brief (not merely finite), and this, plus the lack of a 
stationary phase approximation, requires integration over both fre- 
quency and wave number. These two integrations make near-field 
analysis much more difficult than far-field analysis. 

2. A preliminary near-field analysis is made using coordinates 
fixed in a permanent pulsating source (a harmonic source) . This 
source drifts with the fluid at the center of a uniform circular 
cylindrical jet. For a Mach number of 0.5 and a Strouhal number of 
0.2 the streamwise extent of substantial disturbances is appreciably 
greater than for a Mach number of zero. 

3* The near field is studied for a sequence of transient sources 

at the center of a uniform circular jet. Besults here must be con- 
sidered preliminary since a triple integration is performed numer- 
ically. The integration seems feasible, but the accuracy of the 
methods has not yet been completely explored. The preliminary 
results are reasonable. Most interesting is the appearance in the 
near field of the trend toward Increased downstream disturbances 
which is so characteristic of the far field. 
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Appendix A 


A TWO- DIM EN SIONAL CHAtni^. FLOW SHOWING RESON ANT AND QUASI -RESONANT PEAKS 
Introduction and Discussion 

It is convenient to represent a pulsating source in a jet as the 
sum of components having all possible wavelengths in the streamwise 
direction. Some of these components are true waves, and some decay ex- 
ponentially in the z-dlrection, normal to the flow. Each of these har- 
monic wave forms (infinitely long in the flow direction) moves at its 
own velocity relative to the source. The source itself moves with the 
local fluid, so the velocity of the waveform with respect to the fluid 
in any part of the idealized jet can be determined. Conversely, the 
velocity of the fluid relative to the waveform is known. 

If we then approximate the fluid motion in some part of the jet 
by a uniform velocity flow it is tempting to study for the purpose of 
Illustration a conventional channel flow where one side of the channel 
consists of a prescribed wavy wall. Some boundary condition must now 
be prescribed for the other wall of the channel. It is not possible 
to say precisely what this condition should be. However a generalized 
boundary condition of the form a0 (x, z^) + b0 (xt-y, z. ) = 0 may be 

used. If a = 0 a flat rigid wall exists at z = z^ . If b = 0 
a constant pressure surface exists at z = z^ . Other combinations 
of a, b and y should correspond to the conditions existing in a jet 
though these combinations are not easily found. 

The energy flow across the outer boundary of the channel is ex- 
amined as a function of a non-dimensionalized wavelength in the flow 
direction. Various outer boundary conditions are considered. In some 
cases true resonant peaks appear (the entire denominator passing through 
zero). In other cases quasi-re soneint peaks are evident (the dominant 
part of the denominator passing through zero) . This highly simplified 
model therefore does show features resembling those appearing in the 
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complete shear layer jet, although the peaks appear in a plot of energy 
flow vs, wavelength rather than in mean-square pressure vs. far-field 
angle 0, 

The above results are Interesting and perhaps useful. However it 
must he emphasized that some very important features of the jet, such 
as the jet velocity relative to the ambient air and the velocity of 
the acoustical source, are concealed in the generalized boundary con- 
dition. This procedure could not readily be used for any quantitative 
prediction of the appearance of peaks in the idealized shear-layer 
jets which we are studying, but is possibly of value as a qualitative 
illustration. 

Development 

The development is here outlined very briefly since the problem 
is comparatively simple. Only supersonic flow is considered, this being 
the case of primary interest. 

The origin of the x, z coordinate system is fixed in the mean plane 
of the channel's prescribed wavy wall (see Fig, lA), The flow direction 
is X, the normal direction z. The flow obeys the equation 

_1)C5 _ 0 =0 

''^xx '^zz 

where M is the undisturbed constant Mach number (M > l) and 0 is 
the velocity potential. The prescribed wavy wall is at z = 0, the 
outer boundary of the channel is at z = z^. Let 

= A sln(k^x) at z = 0 (A-2) 



(A-1) 


and let 


a 0^(x,z^) + b 0^(xty,z^) =0 at z = z^ 


(A-3) 


where A prescribes the amplitude and k^ 


the wave number of the wavy 
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wall, a and b are arbitrary constants and y permits an arbitrary 
phase displacement between and at z = 

The following expression for 0 satisfies the partial differential 
equation (A-l) and the boundary conditions (A-2) and (A-3) : 


= — - — ( f ) sin(k^x) sin(k 2 z) + £ sin(k^x) cos(k 2 z) 


abk \ 

-I cos(k. /) cos(k. x) cos(k„z) 1 

^2 / 


where 


^2 = 


* ^ f cos^(k^y) sin^(k2Z^) ! 

[ ]^ = ^(ak^/k^) cosCk^z^) - b sin(k^y) slnCk^z^)]^ (A-5) 

2 

k abk 

[ ]3 = ^ ) sin(2k2Z^) sin(k^y) cos(2k2Z^)]^ 


... 2 / 


(A-4) 


It is convenient to exajulne the energy flow (relative to coordin- 
ates fixed in the fluid). This is independent of z and is given by 
(-yoMc^ jZl) , p being the mass density of the fluid and c the speed 
of sound. This is the local value, and any convenient z (e.g. z = O) 

may be used. If we define by E the mean energy flow per unit area, 
then 


_ A^ab k. ^ cos(k^y) /bMc 

E = 1 ^ 1 

2 k^^ t ! 


(A-6) 


where { ! is defined in Eq, (A-5) and b may be set equal, to unity 
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without loss of generality. 

For a specific example let M = >JZ and =0. Then = k^‘ 

and 


E* = 


2E 

A^/>Mc 


a^ cos^(k^z^) + sin^(k^z^) 


(A 


For a = 1 the flow pattern is precisely that which occurs for an 
infinite fluid flowing past a wavy wall, and E* = 1. If a is large 

p p 

(e.g. a = 10) then a cos (k^z^) is generally the dominant term in 
the denominator, hut passes through zero for (k^z^) = ■^, etc. 

Thus quasl-re sonant peaks appear near these (k^z^) values. If a is 

p 

small (e.g. a = O.l) then sin (k^z^) is generally dominant hut quasi- 
resonant peaks appear near (k^z^) = 0, tr, Z-n, etc. 

As a -♦ oo the outer boundary of the channel approaches a constant 

/ N IT 3tt 

pressure surface and true resonances at “ 2’ 2 ’ 

approached. As a -►O the outer boundary becomes a plane rigid wall, 
and true resonances at (k^z^) = 0, rr, 2 tt, etc., are approached. 



Appendix B 


THE HEATED THREE-LAYER T¥0 -DIMENSIONAL JET 


The problem considered is that sketched in Fig, 13 i but the three 
layers shown now have different temperatures and densities. Outside the 
jet (Region (zj) , the ambient air has temperature T , corresponding 
density p , and speed of sound c . In the outer layer of the jet (Region 
(^) these quantities are T^, p^, c^, and in the central high velocity 
layer (Region (^) they are T^, p^, Cq. 

If the velocity of the jet in Region is and the velocity in 
Region is U^, then the corresponding Mach numbers indicated in Pig. 

13 are defined with respect to ambient temperature; 


"i - V 

”o ■ V 


c 

a 


c 

a 


(B-l) 


Standard thermodynamic relationships give also 



^a 




(B-2) 


We further define 


01 • = Ol/c 

' a 

= k^/w' , = k^/w' 



(B-3) 


If the einalysis described earlier in the report is modified to allow 
for the different temperatures of the three layers, Eq. (ll) becomes 
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|_|2 _ 

where the previous definitions of N, D^, Dg given hy Eq.# (12) for the 
symmetric case and Eq.. (l^) for the anti-symmetric case retain their 
previous forms but the parameters therein defined by Eqt (13) for the 
unheated jet now become 

<r' = "'(ag-a^)^/ 

V = «'(Zg-a^)y 

(B-5) 

(T^/T^)7i:i-K^Mp^-K^ 

(1-M^K^V 

. (TqAi) V(c^/co)^[1+K^(Mq-M )]2-k2 

A*2 = • •• — 

[l+K^(MQ-M^f y c2-(c^/cp2 


and a, defined for the unheated jet in Eq. (8), becomes 

a = a,'a^y(c^/cQ)^[l+K^(MQ-M^)f-K2 (B-6) 

Comparison indicates that the statements made for the unheated jet 

2 

are still applicable. That is, for very large K, is small. 

then is normally the dominant term of the denominator for large K , and 
spikes in the far-field pressure distribution will be found at the zeros 
of . These occur for the symmetric case when 

tanh<r' = ’tana (B-?) 
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which is identical in form to Eg, (l5a), but c* , /* 2 *» ® 

defined as in Eqs. (B-5) i (B-6); the anti-symmetric case leads again 

to Eq. ( 15 b) in form but with the modified parameters above. 

Analytically, then, the heated jet case appears little different 
from the imheated case . The magnitude of I S I ^ is decreased by the 
numerator factor (T /T.) . The location of the spikes (in 0) and the 
range of values of where spikes may be expected has changed (i.e. 

Fig. I4b must be modified to accomodate the new definition of a and 
K > (c /c.) ), but otherwise the problem is qualitatively similar to 

the unheated case. 

When one considers an actual numerical example, however, the shift 
in the K^-range where spikes axe possible can be of major significance. 

The occurrence of spikes for this configuration requires, as pointed 
out earlier in the report, that a be real (corresponding to a wave- 
type solution in the central high-velocity part of the jet) when |k^| 
is large* (exponentiad-type solutions in the outer layer of the jet). 
Heating the jet can have the effect of shifting the K^-range where a 
is real away from the very high values of | K^j . Fewer spikes may then 
be encountered, or possibly only pressure peaks of moderate magnitude. 

An example may serve to clarify these points. Consider a jet with 

Mq = 2.9 a^/a^ = 0.6777 

= 1.5 Sy = 0.2 

For the unheated jet (Tq=T^=T^) a is real for 180° > 0 > 121,76° (and 
also for 104.86° > 0 > 0, but this part of the 0-range produces no spikes 
and so is not of interest here) . Within the range of interest o in- 
creases with increasing 0 from a = 0 at 0 = 121 . 76 ° to a = °o at 
0 = 131 *81° (the Mach angle corresponding to = 1.5 ); « then 

decreases with further 0 increase to a = 228.61° at 0 = 180° . Within 
this a-vaxiation, an infinite number of solutions of Eq. (B-7) Is possible. 


* As before, numerical solutions are considered only for the plane y = 0, 
where K 2 = 0 and K = 1 | , 
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However, if the jet Is heated so that = 3T , = 2T , the 

velocities remaining as atove, it is found that o is real for 180° > 

0 > 148,89° (and for 102.47° > 0 > 0, which is less of Interest), In 
this rainge 2 < < 3*01 *• That is, takes on only moderate values, 

and only moderate pressure peaks are encountered. Also, o increases 
with increasing 0 from o = 0 at 0 = 148,89° to a = 63 . 60 ° at 
0 = 180 °, With this limited range of o -variation permitted, not more 
than one solution (if any) of Eq. (B-7) is probable. 

Thus a velocity distribution which resulted in an infinite number of 
spikes in the fax-field pressure distribution for the source in the \m - 
heated jet has at most one somewhat moderate pressure peak when the jet 
is heated with the temperature ratios indicated. This change is probably 
related to the fact that the actual (local) Mach number of the jet, 
based on local temperature, is greatly reduced. (Of course, one can 
find examples where the Mach number and temperature ratios are such that 
multiple spikes will occur for the heated jet also.) 


* Recall that = cos0/(l+M^ cos0). 
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